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calculated  for  several  model  geometries.  The  calculated  signatures 
were  compared  with  measured  signatures  and  good  agreement  was  obtained 
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CHAPTER  I 
INTRODUCTION 


In  past  years,  research  has  provided  a  great  deal  of  information 
about  the  mechanics  of  scattering  of  an  electromagnetic  field  from 
the  sea  surface.  Much  of  this  work  has  been  motivated  by  the  problem 
of  detection  of  ships  at  sea.  This  report  uses  a  simplified  ship-sea 
model  to  obtain  some  insight  into  the  interaction  between  the  sea  and 
a  ship  when  they  are  illuminated  with  a  radar  signal. 

The  simplified  sea  model  used  in  this  report  is  a  one-dimensional 
composite  model  developed  during  a  previous  study  [1,2],  in  which  the 
method  of  moments  and  physical  optics  were  used  to  calculate  the 
scattering  from  the  surface. 

This  report  expands  on  this  previous  work  by  adding  a  metal  flap 
to  represent  a  simplified  model  for  the  ship.  Since  this  combined 
ship-sea  model  is  too  large  for  the  method  of  moments  approach,  this 
report  further  develops  the  physical  optics  approach  to  study  the 
scattered  fields.  The  first  section  describes  the  surface  and  sets 
up  the  geometry  of  the  model.  Then  a  general  equation  for  bistatic 
scattering  from  the  composite  surface  is  developed.  The  subsequent 
sections  treat  in  detail  the  principal  scattering  mechanisms, 
namely;  scattering  directly  back  from  the  composite  surface;  scat¬ 
tering  from  the  surface-flap  corner;  and  scattering  from  the  edge  of 
the  flap.  In  this  report,  only  horizontal  polarization  is  considered. 


niAPiiR  n 

COMPOS  1 T I  SI  A  SUM  AH  -  SHIP  MOPI  l 


The  purpose  of  this  section  is  to  describe  a  simpl if ied  one- 
dimensional  model  for  a  ship  on  the  ocean.  The  ocean  is  represented 
bv  a  composite  surface  made  of  aluminum  roofing.  A  large  scale 
sinusoid  is  used  to  model  the  gravity  waves,  and  small  trapezoids 
superimposed  on  the  sinusoid  are  used  to  model  the  capillary  waves. 
The  ship  is  modeled  bv  a  flat  metal  plate  placed  at  one  end  ot  the 
model  sea  surface,  isee  figure  1-1 '. 


GRAVITY 


CAPILLARY 


FLAT  PLATE 


Figure  1-1.  Composite  sea  surface-ship  model 


for  the  actual  surface,  designed  to  simulate  the  sea  surface  at 
S-band.  the  sinusoidal  gravity  waves  have  a  poak-to-peak  height  of 
T . b  in.,  and  a  mechanical  wavelength  of  4.1  ft.  This  provides  a  moan 
sguare  slope  of  about  10’,  representing  a  relatively  calm  sea  (Figure 
l-.'l.  The  modeled  capillary  waves  are  double  trapezoids  each  one 
0  cm  long  and  separated  bv  1  cm  (see  F  igure  1-.T1,  recurring  every 
CO  cm  along  the  sinusoids.  This  composite  surface  model  can  be  easily 
represented  mathemat i c a  1 ly .  facilitating  the  theoretical  analysis  of 
the  scattered  fields. 


20.  32  cm 


Figure  1-3.  Geometry  of  double  trapezoid  capillary  wave  profile. 


At  one  end  of  the  surface  model  is  an  eiuht  foot  high  metal  plate 
representing  the  ship.  This  metal  plate  makes  an  angle  with  the 
surface  such  that  a  corner,  not  necessarily  right-angled,  is  formed. 
The  major  contribution  to  the  backseat teri no  from  such  a  configuration 
is  the  double-bounce  reflection  from  this  corner.  Tor  an  actual 
ship  at  sea,  there  are  a  multitude  of  double  and  triple  bounce 
corners  on  the  structure  of  the  ship  itself.  Ivon  so,  the  lamest 
simile  contribution  to  the  return  signature  will  probably  still  be 
the  corner  where  the  ship  meets  the  sea  surface. 

The  ship-sea  model  described  here  has  four  principal  mechanisms 
whi  h  give  significant  contributions  to  the  backseat tered  field. 

These  are  the  backscatter  from  gravity  waves,  backscatter  from  the 
trapezoids,  backscatter  from  the  ship-sea  double  bounce  region,  and 
the  diffracted  field  from  the  plate  edge,  lirst  we  shall  consider  the 
general  physical  optics  formulation  for  this  one-dimensional  geometry, 
and  then  investigate  the  scattering  from  each  of  the  four  mechanisms. 
The  calculated  scattered  fields  will  be  compared  with  the  experimental 
results  obtained  at  the  Encounter  Simulation  Laboratory,  where  the 
ship-sea  model  shown  in  Figure  1-1  was  constructed. 
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CHAPTER  III 

PHYSICAL  OPTICS  EXPRESSION  FOR  SCATTERING 
FROM  A  CYl INDR  I  CAL  SURFACE 


The  first  step  in  obtaining  the  scattered  fields  is  to  calculate 
the  scattering  from  an  arbitrary  cylindrical  surface.  The  approach 
taken  here  is  to  use  the  physical  optics  currents  generated  by  the 
incident  field  and  then  find  the  far  field  radiated  by  these  current 
generators . 


The  general  equations  for  the  physical  optics  currents,  given 
the  incident  field,  were  sunmariced  in  Reference  [1].  The  scattered 
electric  field,  df*s ,  at  a  point  in  the  far-field  of  a  current  element 
Js  is  given  by 


=  jkZ0  (js 

n 


x  rt;) 


x  r'  ds 


(1) 


where  rn  is  the  unit  vector  describing  the  direction  for  R',  the 
distance  from  the  surface  element  ds  to  the  field  point,  k  is  the  wave 
number  and  Z0  =  120”  ohms.  The  surface  current  J  ,  in  the  physical 
optics  approximation,  is 

J.  =  2n  x  ti.  C) 

where  n  is  the  normal  to  the  surface  element  ds  (see  Figure  3-1)  and 
is  given  by 

n  =  -  sin„-  y  +  cos.-  z.  (3) 

The  angle  between  the  normal  n  and  the  vertical  is  defined  to  be 


A 

Y 


Figure  3-1.  Geometry  of  the  normal  to  the  surface. 
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The  magnetic  field  H.  incident  upon  the  surface  can  be  written  as 

-JkRn 

(4) 


fit  ■  'o  (Vh>1/2  ?h  S? 


n 


where  I0  is  a  nominal  current  strength,  and  fv , f h  are  antenna  pattern 
factors,  to  be  discussed  later.  For  an  antenna  of  maximum  gain,  G,, 
which  radiates  power  Py 

!o  =  (PtGt/4„Z0)1/2  .  (5) 


The  polarization  state  of  the  transmitted  field  is  characterized  by  a 
polarization  factor  ?h-  P°r  horizontal  (perpendicular)  polarization 
the  electric  field  £0,  in  the  aperture,  is  in  the  x0  direction 
(see  Figure  3-2).  For  the  geometry  of  the  sea-ship  model,  the  unit 
vector  x0  is  parallel  to  the  coordinate  unit  vector  x.  The  polari¬ 
zation  factor  can  be  expressed  as 


(6) 


Figure  3-2.  Polarization  state  for  the  transmitted  field. 
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mmm 


The  one  way  antenna  power  pattern.  ("(.•  ,n),  is  assumed  to  he  a  product 
pattern  with  one  principal  plane  in  the  x-y  plane  (pattern  angle  st  amt 
the  other  principal  plane  containing  the  i  axis  (pattern  angle  ) 

(see  Figure  3-3),  so  that  F(;'.,n)  can  he  written  as 


=  fy(Ofh(n) 


(7) 


By  assuming  that  f|,(n)  is  Gaussian,  it  will  later  he  possible  to 
reduce  the  two  dimensional  surface  integration  to  a  one  dimensional 
integration. 


2 


X 

Figure  3-3.  Transmitter  and  receiver  coordinate  angles. 
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Combining  Equations  (2)  and  (4)  gives  an  expression  for  the 
physical  optics  currents. 


-jkRn 

Js  =  2  Io  f('-.,n)1/2  -  {n  x  PJ 


Ir 


(8) 


Substituting  the  equations  for  n  (Equation  (3))  and  Pn  (Equation  (6)) 
into  Equation  (8)  gives  11 


Js  -  2I0  f(t.n) 


1/2  e 


-jkRr 


(-sin4>  y  +  cos4>  z)  x 


~(x  x  rj 


x  x  r 


(9) 


From  Figure  3-4  we  have  that 

A  -I  '  A  '  A 

rn  =  R“  (Roro  +  yy  +  hz  +  xx5  oo) 

n 

where 

r  =  sinO  y  -  cosO  z  . 
o 


Substituting  Equation  (10)  back  into  the  equation  for  the  physical 
optics  current  (Equation  (9))  and  simplifying,  we  obtain 


j 


s 


2I0  f(C,n) 


1/2  e 


-jkR 


( Rqcos (0 -<}> )  +  ysin*  -  hcos*t>) 
[(RQcoso-h)2  +  (RQsino+y)2]1/2 


x  . 


(ID 

For  the  antenna  polarization  chosen  here,  the  physical  optics  currents 
are  in  the  x-axis  direction. 

Using  the  expression  (11)  for  3  ,  in  Equation  (1),  and  expandinq 
the  triple  cross  product 
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(x  x  r')  x  r '  =  (  — ,  -  1^  x  +  ---  (y-R ‘ 

n  n  ^r12  J  r,2  o 


cosily  -  — o  (h-R'sim-’)z 

n  i  t  0 

n 


where  the  unit  vector  r'  is  defined  in  Figure  3-4, 

n 


r'  =  —  [(y-R'cos3)y  +  (h-R'sim<)z  +  xx) 
n  D,  w  o  o 

n 


we  obtain 


>  2jkZ 

«es  ■  V2  >0  f(f-n> 


1/2  e 


-jk(Rn  ♦  r;> 


(R0cos(o-<f>)  +  y  si n«i>-h  cos4>) 
[(h-RQcos0)2  +  (Rosino-y)2]1/2 


1  x  +  — _  (y-R'cosH)y  -  — w  (h-R'sini^z 

J  R'  ^  0  R„  0 

n  n 


d  d 

*  JLX 
C0S41 


So  far,  only  the  transmitting  antenna  pattern  has  been  accounted 
for.  The  receiver,  at  the  point  where  the  field  is  being  measured, 
must  also  have  a  power  pattern  and  a  polarization  direction.  We  shall 
assume  in  this  case  that  the  receiving  antenna  has  the  same  pattern 
as  the  transmitting  antenna.  Also,  to  find  the  horizontally  polarized 
component  of  the  scattered  field,  we  assume  that  the  receiving  antenna 
polarization  is  in  the  x-direction.  Equation  (14)  can  be  changed  to 
take  into  account  the  receiving  antenna  by  multiplying  by  the  one  way 
power  pattern,  given  by  (see  Figure  3-3) 
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Figure  3-4.  Incident  and  reflected  ray  geometry. 


fU’ 


,„'),/2  ■  fvU'),/2 


1/2 


(15) 


and  by  taking  the  dot  product  with  the  receiver  polarization 


=  x 


Furthermore,  if  the  maximum  receiving  aperture  of  the  antenna  is 
given  by  Ae,  the  receiver  voltage,  i.e.,  the  voltage  at  the  receiver 
due  to  the  scattered  field  from  ds,  is  found  to  be 


2  j  V 


dVr  "  Tn 


--  f(f.,n)1/2  f(f.\n')1/2  - 


-jk(R 


or 


Rn  R„ 

n  n 


(16) 


* 


(Rqcos(o-4,)  +  ysinj'-hcos*)  d^ 

[{h-R  coso)2  +  (y+R  sino)2]^2  cos^ 


where 


=  y0(Aek2)1/2 


Equation  (16)  gives  an  expression  for  the  horizontally  polarized 
component  of  the  electric  field  at  a  point,  scattered  from  an 
incremental  area  of  the  total  surface.  To  obtain  the  total  field  at 
any  point,  a  double  integration  over  the  entire  surface  is  required. 

Assuming  that  both  the  receiver  and  the  transmitter  have  a 
Gaussian  product  pattern,  then 

1  .  .  2  If,.  2 . 

,  /0  -  o  A  tan  n  -  B  tan  ■. 

f  ( ; . ,  n )  7  2  =  e  2  e  2 

12  12 

,  /0  -  -tr  A  tan  m '  *  k  B  tan  f. 1  ln. 

f(f.',n,)I/z  =  e  2  e  2  ('7) 


where  A  and  B  are  beamwidth  parameters.  The  angles  and  V 

(shown  in  Figure  3-3)  can  be  expressed  analytically  in  terms  of  dot 
products . 


tan;. 


tann 


tanr.' 


rn  *  za  _  y  coso  +  h  si  no 


R.  +  y  si  no  - 

h  coso 

rn  ‘ 

ro 

rn  • 

X 

X 

r  • 
n 

ro 

R  +  y  si  no  - 
0 

h  coso 

rn 

h  coso  -  y  si  no 

Rp  -  y  coso 

-  h  si  no 

(18) 


(19) 


(20) 
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tann ' 


r 1 
n 


R*2  y  cosb  -  h  sim< 
o 


(21) 


Substituting  these  angle  definitions  back  into  Equation  (16)  gives 
the  following  expression  for  the  received  voltage. 


2jVc 
Vr  =  V 


(R  cos(o-s*>)+ysin>,''-hcos^)  /  2 

0  _  _  I  * _ 

y  x  cos4>[  (h-RQcoso  )^+  (y+RQsi  no )  ^  \R^ 


-jk(R;+Rn) 


R 1  R 
n  n 


(22) 


-4x^(--o  *■  -W)  -  !v(ycoso+hsi  no)" /R?1"  -!}(hcos6-ysin6)‘VR?t' 

d  v  d  B,c  t  d  u  c 

(\rj  f\ 

e  d  e  e  d  d 

*  y 


where 


R2  =  Ro  +  y  si  no  -  h  coso 

R'  =  R'  -  y  cosO  -  h  sinO 
2  o  J 


Because  of  the  choice  of  the  two  principal  planes  defining  the  power 
patterns,  we  are  able  to  separate  the  \  and  y  integrations.  By 
integrating  from  minus  infinity  to  infinity  in  the  x  direction,  we  can 
integrate  analytically  the  x  dependent  part  of  the  double  integral  so 
that 


V 


r 


(R0cos(o-^)  +  ysiiH-hcosi ) 

- 

cosvt[ (h-R0coso )  +  (y+RpSino)*"]  ' 


(23) 


* 


D  p  ~ 

-^•(ycoso+hsino)  /R^ 
e  6 


e 


^-(hcosi'-ysintO^/R.^ 


dy 
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where 


_A  /  J_  +  J  2 

-JkW  2U.2  r;2  ' 


.r;2  "  J  V'; 


2  2 


dx 


and 


Rn2  '  R,2  *  <2 


R'2  =  R'2  .  2 

n  1  +  x 


with 


R-|2  =  (h  -  Rqcoso)2  +  (y  +  RQsino)2 


R^ 2  =  (h  -  R^sinn)2  +  (y  -  R^cosb)2 


(24) 


(25) 


Using  the  Fresnel  approximation  in  expanding  the  exponential  phase 
terms 


-jkRn  -jk(R]2+x2)1/2  ,v  'jk(V  2R1  * 
e  =  e  =  e 


-jkR'  -jk(R'2+x2)1/2  'J'k(Rl  +  2RT  } 
%  n  _  i  .v  _  i 

?  -  e  -  e 


(26) 


and  approximating  the  denominator  terms  as 

x_ _ 

2R  2 

a.  e  1 


Rn  (R12+x2)1/2  R1 


3x2 


2R 


,2 


y  e  1 


n3  (Rj2+x*)3/2  R’3 


(27) 
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the  evaluation  of  the  X  integral  in  Equation  (23)  gives 


'  R2R2 


J2" 


-  ( R-j  +R-j 


I  o  T  ?  2  09  0  0  o  o'~o  '  o  o'  "o  o  T  “T*  o  1 /o 

[AR,R^  R^  +AR]  R]  R^+ jkR,  R(  "R^RA  +  jkR  j  R^R  ’  "  +  R( " R;R.^+3RfR$R;"  ] 1 U 


V1  2 


v]  ,'2rv2 


(28) 


Further  evaluation  of  the  integral.  Equation  (23),  would  require  a 
computer.  Notice  also  that  if  Equation  (23)  is  evaluated  by  computer, 
any  desired  vertical  plane  antenna  patterns.  fv(-. )  and  fv(\),  may  be 
used  in  place  of  the  gaussian  patterns  of  Equation  (17). 
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CHAPTER  IV 

BACKSCATTER  FROM  A  MODEL  SEA  SURFACE 


As  pointed  out  in  Section  II,  there  are  four  principal  mechanisms 
contributing  to  the  total  field  scattered  by  the  simplified  ship-sea 
model.  We  consider  first  the  backscatter  from  the  large-scale 
(gravity)  waves. 

A .  Geometry  of  the  Wave  Models 

To  calculate  this  direct  backscatter  from  the  gravity  waves,  we 
first  define  the  surface.  The  gravity  wave  model  lies  in  the  x-ys 
plane.  The  surface  height  is  a  function  of  y  only,  so  that 


h(y  )  =  H  .  +  H  .  *  sin(w  *  y  +  Phase) 
'Js  ei  ei  m  Js 


(29) 


and  the  surface  is  uniform  in  the  x  direction.  The  constant  Hei  is 

the  peak  amplitude  of  the  sinusoid,  wm  =  y-  is  the  mechanical  wave- 

number  of  the  surface.  The  surface  itself  is  40  feet  in  length,  with 
\m  =  4.1  feet  and  He^  =  1.75  inches.  In  implementing  the  compu¬ 
tations  of  the  scattered  field,  the  surface  is  assumed  to  occupy  the 
range  -40<y  -0  with  the  flat  plate  position  at  y  =0. 


The  capillary  waves,  modeled  by  small  trapezoids,  are  superim¬ 
posed  on  the  surface  of  the  gravity  waves.  Each  trapezoid  is  made 
up  of  plane  segments  as  shown  in  Figure  2-3.  The  upper  three  seg¬ 
ments  represent  the  actual  surface.  The  other  two  segments  are 
straight  line  approximations  to  the  original  sinusoid. 


The  scattering  from  such  a  composite  surface  could  be  computed 
from  Equation  (23)  since  the  profile  h(ys)  of  the  surface  is  known. 
However,  in  order  to  simplify  subsequent  calculations,  the  three 
upper  segments  of  each  trapezoid  are  considered  to  be  individual 
scatterers.  That  is,  instead  of  integrating  over  the  surface  of  the 
trapezoid,  we  consider  the  current  density,  Js,  concentrated  at  the 
center  of  each  of  the  segments.  These  current  densities  are  then 
used  to  determine  the  field  scattered  from  the  segments.  If  any 
trapezoids  are  positioned  such  that  a  face  is  in  the  shadow  region  of 
the  incident  field,  then  the  current  density  on  that  segment  is  taken 
to  be  zero.  In  addition,  since  the  field  scattered  from  the  sinusoidal 
component  alone  (i.e.,  no  trapezoids)  is  computed  first,  in  order  to 
illustrate  the  effect  of  the  gravity  waves,  the  total  field  may  be 
found  by  adding  the  field  scattered  by  the  trapezoids  to  this  gravity 
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wave  field.  Thus  the  fields  scattered  by  the  two  lower  segments  of  the 
individual  trapezoid  must  be  subtracted  from  the  total  field  to  account 
for  the  fact  that  the  interval  of  the  sinusoidal  profile,  h(ys), 
occupied  by  the  trapezoid  has  already  been  computed  in  the  gravity  wave 
program. 


B.  Discussion  of  Computer  Code  for  Backscatter 
from  the  Sinusoidal  Surface  (Gravity  WavesJ 

Since  the  profile  of  the  surface  is  known,  the  voltage  at  the 
receiver  due  to  the  field  scattered  from  the  model  surface  is  found 
by  a  numerical  integration  of  Equation  (23).  The  Equation  then  takes 
the  form 

R?RiR 

— 5 — -  (cosO+sinOtanf+y/R  tan4>-h(y  )/R  )  (30) 

0  ‘  S  0 


Vr  = 


-JOc 

2n 


ts2 

1 


si 


-1/2  “3^ ( Rj +R^ ) 


B 

"2 


2 

(tan  4+ tan 


2c 


Ay 


where 


g  =  AR2Rj  2R  '2 2+AR^R1'  2R2+jkR]  R  ‘  2R2R ’ 2+jkR ‘  R2R2R^2+R  j  2R2R ’ 2+3R2R2R; 


,2 


The  variables  in  Equation  (30)  are  shown  in  Figure  4-1,  and  are  defined 
as  follows. 

The  boresight  distance,  R0,  is  the  distance  from  the  radar  to 
the  y-axis,  measured  along  the  radar  axis.  RQ  is  given  by 

Rq  =  zz/sin(^/2  -  o) 

where  zz  is  the  perpendicular  height  of  the  radar  above  the  y-axis. 

The  return  boresight  distance  RA  is  equal  to  R0  since  the  radar  is 
monostatic.  The  incremental  value.  Ay,  is  the  width  (measured  hori¬ 
zontally)  of  the  portion  of  the  surface  contributing  to  the  sum  in 
Equation  (30).  There  is  just  one  physical  optics  current  element  per 
increment.  In  the  program  of  the  Appendix,  y  is  taken  to  be  0.01. 

8y  taking  the  origin  of  the  x,y,z  coordinate  system  to  be  at  the 
intersection  of  the  boresight  with  the  y-axis,  the  distance  between 
the  radar  and  the  surface  point  at  Ay,  R1 ,  can  be  written  as 
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(31) 


R)2  =  +  Ro  sin°)2  +  (h(ys)  -  R0  COSO)2 

=  (y  -  yQ)2  +  (Mys)  -  i>0)" 


where  o  is  the  angle  between  the  vertical  (z)  axis  and  the  boresight, 
and 


y„  =  -  R  si  no 
o  o 


h  =  zz  =  R  coso 
o  o 

The  return  distance  between  the  radar  and  Ay  is  Rj  and  is  equal  to 

Rj  for  direct  backscatter.  The  quantity  R?  is  the  projection  of  R, 

onto  R  and  is  defined  as 
o 

R2  =  Rq  +  y  sino  -  h(ys)  cos 0  .  (32) 

R^  is  the  projection  of  Rj  onto  R^  and  is  given  by 


R£  =  Rq  -  y  cos;'  -  h(y$)  sino  (33) 

where  the  angle  o  is  measured  between  the  horizontal  (y)  axis  and 
the  return  boresight.  The  quantities  R9  and  R9  are  also  equal  for  a 
monostatic  radar. 

The  portion  of  the  surface  contributing  to  the  sum  is  defined  by 
the  limits  y=Ysi  and  y=Ys2-  These  limits  are  determined  by  the  angles 
at  which  the  incident  field  is  20  dB  down  from  the  boresiqht  level. 
Since  we  have  defined  the  point  y=0  to  be  at  the  intersection  of  the 
boresight  with  the  y-axis,  the  value  of  Ysi  will  be  negative  and  YS2 
positive.  The  surface  profile  function  h(ys)  is  always  evaluated 
with  respect  to  the  origin  of  the  x,yc,z  coordinate  system,  that  is 
ys=0  at  the  end  of  the  surface  model  (see  Figure  4-2).  We  then  have 
for  the  limits  of  integration 
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Figure  4-2.  Geometry  defining  the  limits  of  integration. 


Ysl  =  yo  +  ho  •  tan  (0-1-3  *  e3dB) 
Ys2  =  yo  +  ho  ‘  tan  (8+1'3  *  e3dB ^ 


(34) 


where  03^  is  the  3dB  beamwidth  in  degrees,  which  for  our  gaussian 
pattern  nas  been  taken  to  be  15°. 

The  specified  section  of  the  computer  code  in  the  Appendix 
utilizes  Equation  (30)  with  h(ys)  defined  by  Equation  (29)  to  calcu¬ 
late  the  receiver  voltage  due  to  the  backscattered  field  from  the 
sinusoidal  surface  model,  as  the  radar  moves  from  the  end  of  the 
surface  at  ys=-40  to  the  flat  plate  at  ys=0.  This  portion  of  the  code 


requires  input  values  of;  the  height,  D1  in  feet,  of  the  radar  above 
the  y  axis  at  ys=0,  the  look  angle,  «,  of  the  radar  as  measured  from 
the  trajectory  to  boresight  in  degrees,  and  the  dive  angle,  if,  of  the 
radar  in  degrees  from  the  horizontal.  The  output  is  the  voltaqe  at 
the  receiver  (complex  valued)  in  units  of  volts  (see  Figure  4-3). 


C.  Discussion  of  Computer  Code  for  Backscatter 
fronT Trapezoids  (Capillary  Waves)- 

The  most  direct  way  to  include  the  physical  optics  scatter  from 
the  trapezoids  would  be  to  superimpose  their  profile  on  the  profile 
given  by  h(ys)  for  the  gravity  waves,  and  apply  Equation  (30).  How¬ 
ever,  it  is  possible  to  generate  a  set  of  equations  that  adequately 
describe  the  scattering  from  a  single,  arbitrarily  placed  trapezoid, 
so  that  one  set  of  equations  could  be  used  in  a  computer  code  to 
calculate  the  backscatter  from  all  of  the  trapezoids. 


The  field  scattered  from  an  arbitrary  trapezoid  is  found  by 
considering  each  segment  of  the  trapezoid  separately.  The  field  contri¬ 
butions  from  the  three  upper  segments  are  added  together.  Because  the 
field  from  the  entire  sinusoidal  surface  has  already  been  calculated, 
the  scattering  from  those  portions  which  are  now  covered  must  be 
subtracted  from  the  total  field.  For  this  reason  the  field  contribution 
from  the  bottom  two  segments  is  subtracted  from  the  field  ~f  the  upper 
three  segments.  The  net  field  contribution  from  the  five  segments  ;s 
added  to  the  field  scattered  from  the  gravity  wave.  This  process  is 
repeated  for  each  trapezoid  on  the  surface.  When  the  radar  is  moved 
to  a  new  position  a  new  value  for  the  received  complex  voltage  is 
calculated.  In  this  manner,  a  backscattering  S’onature  for  the  model 
sea  surface  is  obtained  for  the  specified  radar  trajectory. 

To  find  the  net  backscattered  field  due  to  a  trapezoid,  we  con¬ 
sider  first  the  field  scattered  from  one  side  of  the  trapezoid  by 
using  a  simplified  form  of  Equation  (23).  Because  each  side  of  the 
trapezoid  profile  is  only  a  small  fraction  of  a  wavelength  long,  we 
can  further  simplify  the  calculation  of  the  scattered  field  by 
considering  the  entire  physical  optics  current  to  be  concentrated  in 
the  center  of  the  segment.  Thus,  no  integration  in  the  y-direction 
along  the  segment  is  needed.  Instead  we  simply  multiply  by  the 
effective  width  of  the  segment.  Also,  since  we  are  dealing  with  a 
monostatic  radar,  the  transmitted  and  received  distances  are  equal, 
eliminating  the  primes  on  variables  which  denote  the  scattered 
distances,  and  doubling  the  transmitted  distances.  Equation  (23)  can 
now  be  written  in  a  form  suitable  for  calculations  of  the  field 
scattered  from  the  trapezoids. 


jV  L(n)sino(n)  (y+R  sino)  +  cosa(n)  (h-R  cos:>) 
v  (n)  =  -  ° - — - 5 - j - A- -  dwy(n) 

coso(n)[(h-RQcostir  +  (y+RQsino)] 

?  ? 

-B(ycoso+hsino)  /R? 

*  e 

-j2kR] 

*  - - p-  e  dx  (35) 

R  ‘ 

where  n=l,**«,5  depending  on  which  segment  of  the  trapezoid  is  being 
considered  (§ee  Figure  4-4);  l(n)  is  either  +1  or  -1,  depending  on 
the  normals  n. 


CENTER  AT  y 


h  ( ys ) 


Figure  4-4.  Profile  of  trapezoid  showing  normals. 


=  L(n)  sino(n)  y  +  coso(n)  z 

L(n)  =  -1  for  n  =  1 ,2,4,5 
L(n)  =  +1  for  n  =  3  . 


(36) 


The  parameter  o(n)  is  also  a  function  of  the  segments  of  the  trapezoid 


Hero,  i  is  t ho  angle  describing  the  slope  of  the  surface  h(ys)  at  the 
midpoint  of  the  trapezoid,  and  is  the  anqle  that  defines  the  slope  of 
the  sides  of  the  trapezoid  (see  I iqure  4-4).  Since  •  > \  detines  the 
slope  of  the  i-th  segment  of  the  trapezoid,  the  width  of  the  i-th 
segment  projected  along  the  y-axis  will  he  a  function  of  , ,  so  that 

dwy(  1  )  -  0.0451. cos ) 
dwy  ( ? )  -■  0.41b-  cos  ( ■  ) 

dw.y(3)  0. 048 -cos (c -i  )  (.18) 

dwy (4)  0. Ob;’- cos ( i  ) 

dwy(b)  -  0. Ob.1- cos ( i  ) 

By  carrying  out  the  integration  over  x  in  equation  (3b)  and 
summing  over  the  five  faces  ot  the  trapezoid,  an  equation  for  the 
received  voltage  contributed  from  a  trapezoid  centered  at  position  y 
can  be  written 


jV 


l.  ( n ) s i n«' ( n )  (y ♦  R  s i  no ) +coso ( n )  ( h  - R  ,cos  e ) 


o 


)  ( M  s  x  r* 

n  n-1  cos. '(n)l (h-Rocoso)‘-+(yfR  sine)'  ] 


0 


-i'jkRj  -|l(ycoso*hsino)'  /R 
*  dwy(n)  e  e 

R,"  (.7R,?  ♦  A/R./  ♦  jk/R,)V'' 


(v) 


(34) 


'1 


'1 


where  the  choice  of  sign  (t)  js  such  that  the  top  segments  (n-l,.\3) 
are  (^)  and  the  bottom  two  segments  (n-4,5)  are  (-).  In  this  equation 
y  is  defined  to  be  the  position  of  the  mid-point  of  the  trapezoid  on 
the  surface  h(ys).  All  other  terms  are  as  defined  for  the  gravity 
wave  calculations. 


The  portion  of  the  computer  code  that  calculates  the  field 
scattered  from  the  trapezoids  is  a  subroutine  in  the  gravity  wave  code. 
This  subroutine  requires  the  position  ot  the  radar,  and  returns  the 
voltage  at  the  receiver  due  to  the  field  scattered  from  the  trapezoid 
superimposed  on  h(ys).  This  complex  voltage  is  added  to  the  voltage 
due  to  scattering  from  the  gravity  waves  so  that  the  program  gives, 
at  each  position  of  the  radar,  the  receiver  voltage  due  to  the  field 
scattered  from  the  composite  surface  model.  The  following  section 
gives  the  results  of  some  typical  calculations. 


0.  Results  for  Hackscatter  from  Surface  Model 

The  signature  for  the  sea-surface  model,  (i.e.,  the  receiver 
voltage  as  a  function  of  the  distance  YV)  as  calculated  using  the 
physical  optics  current  approximation,  is  shown  in  figure  4-5  for  two 
different  trajectory  geometries.  Two  signatures  are  shown  for  each 
configuration,  one  modeling  the  sea  as  a  gravity  wave  alone  (i.e.,  as 
h(ys))  and  the  other  as  the  composite  surface  including  the  trapezoidal 
capillary  waves.  In  all  signatures,  the  magnitude  of  the  receiver 
voltage  is  in  dll  below  reference  voltage  V0,  which  has  been  taken  to 
be  unity.  All  signatures  are  plotted  so  that  they  are  normalized  to 
the  highest  value  in  the  signature. 

Note  that  the  radar  system  used  to  obtain  the  measured  signatures 
includes  a  range  gate  (short  or  long!  to  eliminate  wall  reflections. 

The  effect  of  this  range  gate  has  been  included  in  the  calculated 
signatures  by  inserting  a  gating  function  C(R],1)  in.  e.g.  Equations 
(.10)  and  (.14).  The  argument  1  in  C(R],1)  denotes  that  the  short  range 
gate  was  used. 

figure  4-5a  shows  clearly  the  effects  of  the  periodicity  of  the 
sinusoidal  surface  model .  When  the  contribution  of  the  trapezoids 
is  added  in,  the  signatures  take  on  a  more  random  appearance,  but  the 
underlying  periodicity  is  still  noticeable  (figure  4-5b).  The  lobe 
at  the  end  of  the  signatures  of  figure  4-Sa  is  the  result  of  termin¬ 
ating  the  surface  at  y  =0. 


CHAPTER  V 

DOUBLE-BOUNCE  SCATTER  FROM  A  SHIP  MODEL 


A .  Geometry  of  the  Ship  Model 

The  next  scattering  mechanism  to  be  considered  is  the  double-bounce 
corner  formed  by  the  intersection  of  the  ship  and  the  sea  surface.  In 
the  Encounter  Simulation  Laboratory  experiment  [3],  the  ship  is 
modeled  by  the  8'  high  plate.  To  simulate  different  roll  angles  for 
the  ship  the  plate  angle  u  can  be  varied.  The  radar  approaches  in  the 
yz  plane  from  the  negative  y  direction  along  a  trajectory  defined  by 
the  dive  angle  ^  (see  Figure  5-1). 


Figure  5-1.  Geometry  for  double-bounce  scatter. 
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B .  Simpl if ications  to  Model 

There  are,  in  fact,  two  double-bounce  paths,  one  in  the  direction 
transmi tter-sea  surface-plate-recei ver  and  the  other  in  the  direction 
transmi tter-plate-sea  surface-receiver.  In  both  cases  the  scatter  from 
the  sea  surface  is  in  the  forward  direction.  Because  the  dominant 
scattering  effect  from  the  sea  surface  is  the  forward  scatter  from  the 
large  scale  structure,  h(y-),  we  can  neglect  the  effects  of  the 
capillary  waves  in  calculating  the  double-bounce  effect.  This  will 
greatly  simplify  the  field  calculations  and  the  computer  code. 

Even  after  neglecting  the  capillary  waves,  a  direct  calculation 
of  the  scattered  fields  would  take  an  unreasonable  amount  of  com¬ 
puter  time.  Instead  we  may  take  advantage  of  certain  symmetries  of 
the  problem  to  further  simplify  the  calculations. 

For  any  configuration  of  the  sea-pl ate-transcei ver  positions,  the 
geometry  can  be  made  symmetric  about  the  flat  plate  by  using  image 
theory.  That  is,  since  the  plate  is  an  ideal  reflector,  we  may  image 
the  surface  and  transceiver.  By  appropriately  defining  the  limits  of 
integration,  we  can  then  effectively  remove  the  plate  (see  Figure  5-2). 
This  can  be  done,  of  course,  only  if  the  sea  surface  behind  the  plate 
can  be  neglected.  For  the  study  of  the  double-bounce  scatter  this  is 
indeed  the  case.  Reflection  off  the  plate  is  taken  into  account  by 
multiplying  the  field  passing  through  the  plane  of  the  plate  by  a 
minus  sign  (for  horizontal  polarization). 

Another  simplification  is  that  the  two  ray  paths  from  transmitter 
to  receiver  through  any  point  C  on  the  plate  are  of  egual  length. 
Therefore,  by  reciprocity,  we  may  calculate  the  field  due  to  a  first 
reflection  off  the  sea,  and  multiply  it  by  two  to  account  for  the 
ray  hitting  the  plate  first.  The  resulting  confiouration  is  shown 
in  Figure  5-3. 

Thus,  by  image  theory  and  reciprocity,  the  problem  has  been 
reduced  to  a  simple  application  of  the  equation  derived  in  Section 
III  for  bistatic  scattering.  For  a  horizontally  polarized  field 

2jV  RR' 

Vr  =  -  -  -j—  (R  cosn+R  sinntanj>+ytan^-h(y  ) ) 

2n  \  1 
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Figure  5-2.  Geometry  resulting  from  image  theory. 


where 


g=ARjRj  ‘R^+AR^Rj  2R2+R^R2R*  "+3Rj  R^SjkR^R^Rj  R^+R]  R ' c ) 
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4.2. 
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and  h(ys)  is  the  sinusoid  defined  by  Equation  (241.  The  appropriate 
limits  of  integration  will  be  discussed  in  the  next  section. 

C .  Discussion  of  Coniputer  Code  for 
Double-bounce  Scatter 


The  field  which  is  reflected  back  to  the  radar  by  the  presence 
of  the  plate  can  be  calculated  by  a  numerical  integration  of  Equation 
(40).  By  reducing  the  integral  to  a  finite  summation  over  small 
increments.  Ay,  we  have 


2JV0  Ys2  R2R' 

Vr  = -  }  (R0costHRosinetan,*.+ytan4.-h(y  )) 

2.  Ys1  1 


1/?  -jk(R,+R')  -  |  (tan^.+tanS’. 1 ) 

*  g  '  e  1  e  £  Ay  (41) 


where  g  is  defined  in  Equation  (40).  Unlike  the  formulation  of 
Equation  (30)  for  the  scattering  due  to  the  gravity  waves  alone,  the 
printed  variables  and  the  unprimed  variables  are  not  equal.  This  is 
because  of  the  simplifications  introduced  by  the  image  theory.  Figure 
5-3  defines  the  variables  for  this  scattering  equation.  The  terms 
R]  and  Rj  are  the  distances  from  the  incremental  segment.  Ay,  on  the 
surface  h(ys)  to  the  transmitter  and  receiver,  respectively.  These 
distances  are  given  by 

Rl2  =  (*-V2  +  (h‘ho)2 

2  f  (y.y+y0)sin*  \2 

=  yy-y  +  - - - +  (h(y_)  -  Rising  (42) 

1  V  sin(J  -2u+o)  J  s  0 


where  the  angles  and  the  quantities  y,y0,yy  and  h0  are  as  shown 

in  Figure  5-3.  The  receiver  boresight  distance  R^  is 


(43) 


R ’  =  R  +  (yy+y  )sin2u/sin("/2-2..  +  ’1) 

0  o  •  0 

where  R0  is  the  boresight  distance  to  the  transmitter.  The  antenna 
pattern  angles  •_  and  '  are  as  defined  in  Section  III.  The 
projections  and  RJ,  are  given  by 


R2  =  Ro  +  ys’n  "  h(ys)cos!’ 


R2  =  Ro  '  (y~2y0  "  2yy^  C0S;'  "  h(ys}  s1m'  •  (44) 


The  limits  of  this  numerical  integration,  Ys]  and  Ys?  are  determined 
by  the  20  dS  points  of  the  antenna  pattern,  as  for  the  gravity  wave 
integration,  with  one  additional  constraint.  The  integration  must 
take  into  account  only  those  ray  paths  for  which  either  the  trans¬ 
mitted  ray  or  the  reflected  ray  passes  through  the  plate.  The  limits 
of  integration  must  satisfy  both  conditions  simultaneously  (see 
Figure  5-4). 

The  calculations  of  the  double-bounce  scatter  are  made  in  the 
subroutine  DBS.  The  subroutine  requires  as  input  the  position  of  the 
radar,  and  size  and  position  of  the  plate.  The  subroutine  returns  to 
the  main  program  the  complex  voltage  at  the  receiver  due  to  the 
reflections  from  the  ship-sea  corner  back  to  the  radar.  The  following 
pages  show  typical  signatures  of  several  double-bounce  geometries. 

D.  Resul ts  for  Double-bounce  Scatter 
from  the  Ship  Model 

The  signatures  for  the  double-bounce  scatter  contribution,  as 
calculated  using  Equation  (41),  are  shown  in  Figure  5-5  for  several 
different  model  geometries.  The  magnitude  of  the  receiver  voltage, 
in  dB  below  the  reference  voltage  Vo>  is  plotted  as  a  function  of  the 
position  yy  of  the  radar.  As  before,  a  range  gate  function 
C[(R-|+R^  )/2,l  ]  must  be  included  in  Equation  (41). 

The  effect  of  the  depression  angle  »  can  be  seen  by  comparing  the 
signatures  of  Figure  5-5a  or  5- 5b .  Note  that  for  near-grazing 
incidence  angles  the  double-bounce  contribution  extends  over  a  greater 
range.  The  effect  of  the  dive  angle  „  can  be  seen  by  comparing  the 
signature  with  ,»=60"  and  u  =  5°  in  Figure  5-5b  with  the  signature  in 
Figure  5-5c  that  has  the  same  a  and  u.  A  larger  dive  angle  results  in 
a  signature  less  extended  in  range  and  shifted  in  position  towards  yv=0. 
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Comparing  the  signatures  of  Figure  5-5c  illustrates  the  effect  of  the 
plate  angle  u.  The  radar  return  becomes  stronger  over  a  wider  range 
for  obtuse  corner  angles  (v^O.),  and  decreases  in  magnitude  with 
compressed  range  for  acute  corner  angles  (u^O.). 


Figure  5-4.  Geometry  for  defining  the  surface  of  integration. 
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CHAPTER  V! 

BACKSCATTER  FROM  THE  PLATE  EDGE 


In  the  previous  section  it  was  found  that  a  significant  return 
can  be  expected  from  the  double-bounce  mechanism,  there  may  also  be 
a  contribution  due  to  diffraction  from  the  edge  of  the  plate.  When 
an  incident  or  reflected  field  interacts  with  the  edge  of  the  plate, 
there  is  a  field  diffracted  in  all  directions.  Because  of  the  method 
of  solving  for  the  double-bounce  scattering,  the  edge-sea  surface 
interactions  have  been  accounted  for.  Hence,  the  only  field  contri¬ 
bution  that  has  not  been  considered  is  that  which  is  diffracted 
directly  back  to  the  radar.  The  Geometrical  Theory  of  Diffraction 
(GTD)  is  used  to  study  this  contribution. 


A.  Geometry 

The  GTD  solution  requires  only  a  knowledge  of  the  plate  and  radar 
positions;  the  sea  surface  does  not  affect  the  direct  diffraction  from 
the  plate  edge.  As  before,  the  radar  approaches  from  the  neoative 
y-direction,  with  the  distance  from  the  radar  to  the  plate  edge  defined 
as  Rg.  The  geometry  of  the  plate  and  radar  is  shown  in  Figure  6-1. 

B .  Equation  for  Backscatter 


Since  here  we  are  only  concerned  with  the  edge  diffraction,  the 
flat  plate  can  be  considered  as  the  limiting  case  of  the  wedge,  for 
which  much  work  has  already  been  done.  In  reference  [4]  an  expression 
is  given  for  the  diffracted  field  from  a  wedge.  For  receiver  distances 
R  from  the  edge  which  are  approximately  equal  to  the  distance  R0 
of  the  source  from  the  edge,  this  diffracted  field  is  given  as 


- jk(R+R 


v  i  i  i 

i .  d  "  «|E 


(45) 


where  the  choice  of  sign  (-)  is  determined  by  the  polarization  of  the 
field.  The  minus  sign  applies  for  an  electric  field  polarized 
parallel  to  the  edge,  and  the  plus  siqn  for  perpendicular  to  the  edqe. 
Vb  is  the  diffraction  function  and  is  defined  in  Reference  [4],  The 
angles  g- '  and  ij'g  are  defined  in  Figure  6-1,  and  are  equal  for  a  mono¬ 
static  radar.  Similarly,  the  distances  R  and  R0  are  equal  and  are 
defined  in  Figure  6-1  to  be  the  distance  Rq.  Equation  (45)  then 
simplifies  to 
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field.  The  calculation  is  made  for  each  new  position  of  the  radar. 
Some  diffraction  signatures  are  shown  in  the  following  section. 


D.  Results  for  Edge  Diffraction  From  Plate 

The  signatures  for  the  field  diffracted  from  the  plate  edge  are 
shown  in  Figure  6-2.  The  model  geometries  that  are  shown  for  the 
diffraction  mechanism  are  some  of  the  same  geometries  shown  for  the 
double-bounce  mechanism,  Figure  5-5.  Again,  the  range  gate  function 
C(Rg,l)  must  be  included  in  Equation  (46). 

Note  that  as  in  the  double-bounce  backscat.ter  (Figure  5-5),  the 
diffraction  signatures  extend  over  a  greater  range  but  are  smaller  in 
magnitude  for  near  grazing  incidence  angles.  Tilting  the  plate  results 
in  a  shift  of  the  position  of  the  maximum  return.  Increasing  the  dive 
angle  also  results  in  a  signature  extended  over  a  greater  range. 


CHAPTER  VII 

COMPLETE  MODEL  FOR  SHIP-SEA  SCATTERING 


The  previous  sections  illustrated  the  four  principal  mechanisms 
for  backscattering.  Each  mechanism  was  considered  independently  of 
the  others  to  show  how  each  contributed  to  the  total  ship-sea  back- 
scattering  problem.  Since  we  are  measuring  the  voltage  at  the 
receiver,  superposition  allows  us  to  sum  the  voltages  due  to  each  of 
the  independent  scattering  mechanisms.  The  result  is  the  measured 
receiver  voltage  due  to  scattering  from  the  complete  ship-sea  model 
(Figure  1-1 ). 

A .  Comparison  with  Measured  Results 

Figure  7-1  shows  the  comparison  between  a  signature  calculated 
by  superposition,  and  the  signature  for  the  corresponding  geometry  as 
measured  at  the  Encounter  Simulation  Laboratory.  A  check  of  the 
accuracy  of  the  absolute  level  of  the  calculated  receiver  voltage 
can  be  obtained  by  using  the  return  from  a  calibration  sphere  as  a 
reference.  The  sphere  return  was  calculated  using  the  radar  range 
Equation  [5]  for  a  14"  diameter  sphere  placed  18' 7"  from  the  trajectory 
of  the  radar.  From  Figure  7-1  it  is  clearly  seen  that  the  calculated 
and  measured  magnitudes  of  the  scattering  signatures  are  at  the  same 
level  relative  to  the  calibration  sphere  return.  In  all  of  the 
signatures  in  Figure  7-2,  the  level  of  thp  experimental  signature 
relative  to  the  calculated  signature  was  determined  by  super¬ 
imposing  the  calculated  and  measured  sphere  responses.  The  measured 
signatures  are  also  identified  by  the  experimental  run  number  [3]. 

Figure  7-2  shows  the  comparison  between  the  complete  model 
signatures  and  the  experimental  data  for  a  variety  of  model  geometries. 
In  all  signatures  the  short  range  gate  was  used.  Note  that  in  most 
of  the  signatures  the  double-bounce  contribution  is  easily  identifiable 
as  a  change  in  the  magnitude  of  the  signature.  In  the  geometries  of 
Figures  7-2a  through  7-2c,  the  contribution  from  edge  diffraction  can 
be  seen  as  the  last  lobe  on  the  signature.  In  all  of  the  signatures 
the  quasi-random  return  from  the  trapezoids  is  modulated  by  the  gravity 
wave  (sinusoidal)  contribution  in  the  "clutter"  region  of  the 
signature. 
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CHAPTER  VIII 
SUMMARY  AND  CONCI  US  IONS 


The  purpose  of  this  report  was  to  develop  a  method  for  numerically 
calculating  the  backscatter  from  a  simple  one-dimensional  sea-ship 
model.  The  surface  was  a  "composite"  model  of  the  sea  surface  with  the 
large  scale  gravity  waves  modeled  by  a  sinusoid  and  the  capillary  waves 
modeled  by  trapezoids  periodically  superimposed  on  the  sinusoidal 
surface.  The  ship  model  was  composed  of  a  metal  flap  mounted  on  the 
surface  model . 

For  this  ship-sea  model  there  are  four  principal  scattering 
mechanisms  that  contribute  to  the  backscattered  signature;  namely, 
backscatter  from  the  gravity  waves,  backscatter  from  the  capillary 
waves,  double-bounce  scatter  from  the  ship-sea  corner,  and  dif¬ 
fraction  from  the  edge  of  the  plate.  Since  the  surface  model  was  too 
large  to  employ  the  more  exact  moment  method,  the  first  three  of  these 
mechanisms  were  analyzed  using  the  physical  optics  approximation, 
while  the  last  mechanism  was  calculated  by  a  direct  application  of  the 
Geometrical  Theory  of  Diffraction.  The  voltage  received  at  the  antenna 
terminals  was  calculated  for  each  of  the  scattering  mechanisms. 

The  results  were  presented  in  a  form  showing  the  contribution 
from  each  of  the  individual  scattering  mechanisms,  and  these  were  then 
combined  to  produce  the  total  received  voltage  as  a  function  of  radar 
position.  The  resulting  signatures  compared  favorably  with  the 
measured  data  from  the  Naval  Weapons  Center,  both  in  absolute  level 
and  in  the  spatial  variation  of  the  signature. 

Thus,  the  technique  presented  here  provides  an  accurate  and 
efficient  method  of  determining  the  backscatter  from  simple  one 
dimensional  ship-sea  models.  The  accuracy  of  the  results  were  limited 
mainly  by  the  use  of  the  physical  optics  approximations .  The  runnino 
time  of  the  one-dimensional  surface  programs  could  be  reduced  by 
employing  a  faster  integration  scheme,  since  no  special  effort  was 
made  to  produce  an  efficient  computer  code. 

It  should  also  be  noted  that  the  physical  optics  technique  can  he 
applied  to  a  wide  variety  of  surface  structures  and  experimental  con¬ 
figurations  provided  the  integration  increments  are  properly  chosen. 

For  example,  the  antenna  power  pattern  or  polarization  of  the  radar 
system  can  be  easily  changed  from  those  used  in  this  study.  The 
ability  to  separate  the  contributions  of  the  four  scattering  median 
isms,  and  to  associate  a  range  and  doppler  shift  with  each  point,  on 
the  surface,  provides  a  convenient  method  for  investigating  the  effects 
of  a  variety  of  range  gates,  doppler  filters,  polarizations,  etc.  on 
system  performance. 


APPENDIX 

COMPUTER  CODE  SSSCAT 


A  listing  of  the  computer  code  discussed  in  the  text  is  presented 
here.  To  facilitate  the  understanding  of  the  code,  the  symbols  used 
in  the  program  have  been  used  in  the  text  whenever  possible.  Note  that 
the  code  may  employ  programming  features  not  universal  to  all  computer 
systems. 
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rtul  a  C  t  X  P  (  P  f  P  I 

>1  a  hi'«h?»«>  .  aWP?  «PtPl /PFN1 

I'M  a  A  T  AM  w>  *hF  1*C3S<  »P*Y5*P6iAol  «Pl/l  1*0.  )  ) 

Ajl  a  (' Tr  +  S  I  Ha  TAN  I  pm  )  *Y /Po*  T  AN  i  ph  ) -m/ho 

.'nr  a  A  p  1  /  k  1 

rL*- 5  =  I  »►'«-«/?. »T'f 

rtL  OP  a  t  X*  I -M/2.  »T*ltP?  > 

ftEOP  a  LU  f 

i’  a  (MlapPl >/?, 

i’X  a  -  l  .•Xl*dY«AHt»Cl.LP»rLtOP*YO»C  ( ■».  *t /SOKT  (PT  «i>  .  > 

*  xY  a  V X Y ♦ V  x 
Y  a  Y»i-Y 

IM Y.tf .  YSi  »  GO  IP  31 
”0  91  Jjai.fO 
"O  91  » a  1  ,/ 


T.  i T  S 
T  ■  • 


n:\.  rioABUK. 
.UJU  DO  sot 


lbo  V,-,s-JJ«0.bbab7*0.1;>3 

167  TF(K.tC.;;iY'*aY*-0. 19587 

Ibb  <‘ALU  SPHS(VStFiTMJTWi') 

lb9  VXYsVXt ♦  V  STP 

170  81  ru'.TIiVE 

171  31  Call  lUSliTOt  vrPS.YY.ZZ> 

l7a  wxYY(^>  =  vX»*VufS 

i/3  «i  continue 

179  WK  t  Tt!  (  1  I  (VXYY(N).*'al.KKK  I 

i?3  FALL  GETLP(NCP> 

1 7b  CP  a  ('‘iCP-NCoi/lUfl. 

1/7  *HITe<  J,t.OO>CP 

17b  800  MJHMAT < »u. ,10X« • TOTAL  CP"a«,F8.;>) 

179  CALL  L*  I  Y 

180  r.Yn 

181  C 

i8z  r 

18J  OLCCK  °A  f  A 

18b  I  0OICAL  lUST.I  OEdL'G 

13b  r  1  AENSIO'V  Plb).L<3> 

18b  °t_AL  L 

187  rO  -,«IO‘4/Tt:sT/LTESt  .LCE°UG 

1L8  ru^MO  J/PtS/Pl.TPI.-lPK 

ld9  CO^iOiM/Pf  S/P.L 

190  'ATA  PI.  TPl.(,PP/3. 191926b,  6.  2831 S31 . J ?. 29577950/ 

191  hata  lies! ,Lrjr«uG/. false. ..fal^f./ 

192  I'ATA  P/0. li«Ufe9,0.f>*6itFO,O.(l«'>49,O.l)?f  .0.028/ 

19J  PA  T  A  L/-1.C.-1.U.1.0.-1.0.-1.0/ 

19>»  .  .91 

19b  r 

19b  C 

19  7  r 

198  r  PA/I6E  GATE. 

199  rU-iCTlPM  C(9,“) 

200  C 

«J1  r 

Z'V  r 
203  C 
209  r 
23*>  r 
2 Ob  c 
207  C 
20o  C 
209  r 

2l('  <* 

211  u  Cl  'T 1  iy"E 

212  WETul., 

213  Cm' 

219  c 

21b  C 

21b  r  SUPKOUTlfit  SP»S  CAlCULA/PS  The  SCATTF*1'S 

217  r  IhOM  IMF  SPALL  FP-'oEZOIOf  L  PY  9  I  uP8A  Y  I Otlb 

21b  r  Tuts  19  A  pSAGP  SCaTTkr  ryfil » loUT  ION 

219  r 


sunKouTt.Mt  SP'tS«vsrP.T»4,Tiy/') 

p(5i.Ltst,H<us),'rPiai.r*.TiM.si3(M 

rui*SOii/Sui<r /PNASl  .  m£  1 1 .  mP,po,*m 

,  mo.p.kp 

r  0*M0N/P  t  S/P  l  ,  TP1  ,;'PW 

ro“M  /  yy  *i?  *  y  j.ap 

rO'»«ON/PTS/P.L 

CO^Pct-X  o.VsTP.OC  w.urM  .VbiAPtitCANR 
»tAL  CiUP.Kv.1 
0»t  0. >1.0 ( 

.ajabj/i.ipA;?i 

«TM  =  SINttHU 
I  IM  I  CObt  Tnl  I 
1  isY* 

"l'  =  "fill 

!H  YS.GE  .'SMS  Inf  I*n. 

‘■L.  =  ci 

’'Ms** *Mr  I  i*s  (  w*b  Y  S Sc*P y  / 1  *j|j  #  | 

O'tsAT.Ulli'Ml 


in‘  'OLU'ii'c.  pti  t'ts  fur  PfP“LTnt  pr  Tt-«- 
Tn \P£C'Ii  S.  MS  AMI.  YP  Of F  t AT  Y».|  POSITION  Or  Y*t 
rij’Af  |T  U.lvENT  IM  fAO  StG**rNT. 


iM 

<?a<? 

;aj 

r 

M- 

’37 

i3f' 

4»i'  r 

«bl 

^b.^ 

1** 

«:»«-  r 

<ib7 
^bi 
ibM 
i7u 
it  l 
iTj  r 
^T.x  r 
4  74  r 
.7  S  n 


“Slllsi'll  »«SUl«M.'-fPi> 

••SUIaPI.'  ••CtStfPS) 

“it  ATsPt  A  t  NU*£.MS  1 

••bt'*lak't-)«S|M-tP'!) 

•  it  a  I  s>'  (MbSlMt  »f  ms  ) 

YlM  1  la-Mt  1  )»COSt7Mj»£PS> 

YPt rfl=-MI , )»StMtf MSI 
Y»M  A)  al’l  A)  •(  PS  •  i'U*f  PS  I 
YMtHlaPpObCPStEPS) 
rp( 5|a-MtS)*C0S»LPSI 

c  I »•  <  1 »  :Mhc«i  PS 
S  1»  t  .>  I  at  MS 
SUM  AiaMMi  -I  PS 
si'-.  (•»)=»:  ps 

>  i‘- 1  a  i  sci  s 

*«Y  <  1  )  3.1  Hi- »COSt  RMpbf  PS  I 
'.rtMis.i.ujvapoattPb) 

-'■Yt  '»*.C»P*C0StPH'-tpSl 

:  «Yt«*i  =  ,i  ay«cnsttMS) 

''«Yts|s.i.;.»cnS(CPSI 

Ym's  SFC 1 1 1 CU.COi.ATrS  Thl  VPlTAGC  AT  tmc 
••'CCCTVPP  ICC  TP  Th£  CjltP  ScaTIE«E0  TK(iv 
‘  ACh  ;il  nr  1  Hf  104Prfi>jn.  ivOCTSt 


{IBIS  rWl'VNT  yr;  cr  •• 

V*  .VV|-\  ... 

S t  1  ■  m i c  1  .\ ,  ••-• 

r  “r?pr'!  ’  ■  i  .  v 


in-  rKAcricABi*,, 

a  ■  Tw;.  A 

.'it  r*o  JWT 


276  c 
in 

if  6 
279 

28  o 
281 
2  82 
28  J 
284 
2  83 
286 

287 

288 
2«'-> 

29, 

291 

292 

293 

294 

295 
2*6 
297 
29b 

299 

300 
3<M 

302  7o 

30.  « 

3J4 

303  r 
30b  f 
307  r 

30  >♦  c 

309  r 

si o  r 

311  r 

312  f 

313  r 

310  C 
315  C 
31b 
317 
310 

31  .> 

Hi, 

32 1 

322 

323 

32  4 
325 

324 
327  r 
32<>  r 
32  t  r- 
33o  r 


«STt»*(  9,,0.  I 

ou  73  'isl.S 
uh:HS  |M*k 
Y  T  =  T3*YY*Yh4YP  (  N  I 
F 2»rt0 ♦ Y  T» ST k-mH*C  Tm 
qH2s82»«2 

OKl*tYT-Yp)»«2*«H9.My >».? 

O 1*S0KT ( b  K 1 ) 

T,^s(Tt«crH*MH«SThl  /92 
TiYTs  7 

IMT'Jt.Ot.i.lTMT*?. 

<:=’./  (Fl»A/«K2 
TikP/P 1 

't,,2=C“OLK(S,Tt 
•'tf  IsCSOPTCf  M?) 

0£C5U  <  41  *SI  '  ISIot  4»  )•  1  YT»KC«SrM)*C0S151.",(  jll  *(Hh-«0«cTh» 

‘>tL'S(,tM'bo»CTH>««24<  YT*4o*STul««? 

RLUE=L0S(SIbtft>  »*SjPT(l>L'') 
ri_2JsEXP»-H»T>:T) 

AHl's-2  •  ‘o**-  l 
CA'»t*  =  C£XP(  «►»  I 

05  =  -VO*J»>Jl  t  /PUUl  »0WY(M*ClFP»Ca»K/CFM/2./K81/S2«|  (FI) 
vi=ws»r(4i,'  ) 

IF  (  N.jf.  .4  )  VS  =  -V/S 

F!»TPavbTF*VS 

rof,TM:iE 

PtTUKN 


Tnls  StMrtOL Mr  r  r<H  CllATcs  Tuc  POUttl  F  80.nCt 
c  C  A  T  t  L  9  t-40"  7  FLAT  Pl.ATr  AT  A  SPECIF  I  £  f  a  Y(*(  (.  TO  A 
OJiPOSTTL  SI'iiSJIJflL  SI'RfaLE.  T^t  SOLUTION  ALSO 
ru6T.iriS  TME  FHGE  OIFFPAfTIOP  FROM  Tbt  TIP  OF 
T HE  PLATT.  TwF  FIELH  IS  HORT2UNTALLY  POLaRIZFO. 


PU'iROuT  tut  '0S'»TO  t  Toll’s*  vV  .271 

F’J 990.4/  a1'  6/  T“FTA*  T>(I  .PHlr  «PL!  •  Ai  FHA.  ALP •"  J.F  L9 

ru.'IRON/SuRF  /PlASF  t  wt  1 1  .  mp  ,  L)p  ,  K9 

ro^sof./AYi  /P4v«8»4,a.p,cl.v9,',.kp 

rU"PLt\  l  EA.O.rF'1!  .PEP 

COUPLET  PtPl .  xt « v»  .  VXY  .  vrifS 

C OF 90*4/ 9 1 5/ PI  .7^1  tOPP 

C0MPLlT  EXYY.Un 

real  A'V'Uin'R 

HLt  =  HE  II 

•>TA  =  r,!/2.«TMl*2.«,*U9 

TmT  F  OLLOl  I  NO  SITS  UP  THr  (Apl'NuAR  I ES  OF  INTToRaTJON 
AS  uEFT'itU  PY  Tml  9tA»L!nTM. 
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331 

332 

333 
339 

353 

336 

337 

338 

339 

390 

391 
39/ 

393 

39*+ 

3  «  5 
39b 

397 

398 

399 
330 
351 

35*  r 

354  r 

359  r 

336  r 
35o  r 

337  C 
35., 

35" 

360 

361 

36/  31 

363 

36*+ 

363 
366 
36  7 
366 
36  7 

370 

371 

m 

3  73 
379 
375 
3  76 
577 
3  70 
379 
350 

381 

382 
30  * 

389 

385 


T.+l'S  a  -MPaYY/ (  22-hP  > 

YU  a  -«0»SI6t THI I 

»U  *  22/5171  (PI/?. -TM1  ) 

r>J  a  Pl/*.-2.«i'*urUrHL 

AKP  I  <YY*Y''|«SIN(2.»9li8|/SI*'<dYl 

KKO  a  HO*AKP 

■+K0  a  KIU** 

Pl\Pf)  a  9)*0*»2 
CTM  a  CnsilTn!) 

STH  s  Sir.  (Toll 
MU  a  22 
Yu  -  -M0»STM 
VPUS  a  YT  +  Y'* 

Y31  a  yo*ml  *TAN<  tMl-l.3«o*y//npi() 

“1  a  T‘»l*i  .3*AWV/?0H 

IP1P1.6T.1.«?)«*1S1.55 

v02  a  Yn  +  M(i*TAN(Pl  ) 

v3l9aYP08*Y>  EG 

vG-HaYPOS 

TF«Y91*i.ur.YSlt»51sTSlP 
I F  <  YS/.r,T.YS2M|I32sYS2P 


r-irs  Ol-CTIu11.  INTEGRATES  (HER  Tmf  ANTe';N7  footprint. 
rnr  at9uLT  IS  YH£  vOLTAGF  UFOElvFO  PAC*  aT  T-E 
Pul  MT  IF  Thf.vS»ISSlv)N.  (VOLTS) 

<  xy  a  lo.u.u.ni 
»►  (YSl.PT.Yi.2lG0  to  31 
CY  a  *'.01 
Y  =  TGI 

Yjj  a  .YO-YYYY 
-e.1  S  "Fll 

TE  (YS.*;t.-.r6333M5  tap. 

+•  a  MtI*MtI«SI'"(9,'«YS+PMfit,Pl/jBp,  , 

r+  =  YS+.np- 

•  <*  a  MtI*MtI*sI'l(.'«»YK*PMASr*Pl/180. 1 
I  K'I*Y»STH.H»(.Tm 
R93  a  P2»»2 

a  K 1  S  <  Y-Yt  1  •  •2'M  M-hu)  *«3 
r 1  a  sc« t IPHl) 
r AT  a  YPU£/GI-.(3T  )  •SI>,(PTA  ) 

R*<H1  a  (  Y  Y*0  AT-»  *  •  »2*  (  M-MPO»a  I  *+ (  PTA  I  )  ••a 
<*V!  a  SJKT(P«Pl> 

T.VF  a  (Y*CT»+H«STm)/k; 

Ti+E  2  a  T*  E**2 

*<P2  a  PH  -<  -2,«YY-2.*V('*Y)*roSt8TA)-H»SI  ,(PlA  ) 

9962  a  9+a#«2 

Tuc  P  a  H.CGS(pTAl/RP2»(»2.*YV»^.»Yu  +  Y>*SI’j(PTA)/PP2 
TF«  7  ‘EP.t'T.l  .M  TNEPal.5 
Tt+E P*  a  Txf.P««3 
T  f  (  T  it.P2.uT  .2.1  Tr.-tPaaj  , 

T  a  -XP»  (  91  •°«Pl«882*Ht,P5*'<fi?,KHP2»RP  1  *»ul  ) 
f  a  -<  A»l(8l»6RPl»-i9P2*683*K8Pl  »KBP.>  +  J.aWRl  •PppaMHu*  ) 
S  a  S-A^hXl'PKPl^XRZ 


3»fe  i>tN2  -  C*PL> IS* f ' 

367  i  tMl  =  CSyKl(nC92l 

386  eL  s  -kp» <H 1 *RP1 ) ♦FI//. 

389  RtO  a  0.0 

39u  Pt°  =  C^LX ( K£0.«t I 

391  »tpl  =  CtXFMPEP' 

3^2  >1  -  Ko*92«c  .*Rpi»ptpl/Of  N1 

393  Ph  a  ATrti»(wv«HEI«CrS(*>'*rS)*PMASF.Pl/l80. )  | 

399  P.bl  a  CTh9SlM*T*MP»«l9Y/R0*TAM|pM»-H/K0 

395  7«E  a  AHl/hl 

J9fe  Ti.ro  a  FXP<-B/2.*l9f2> 

397  CLfOP  a  LXP(-B/2.*T9l P?) 

39*  P  a  («1*6P1I/?. 

399  V*  *  *.**l»l'T*At3t*CLt.r«Cl.E0P*VO»C(K.*"»/SQRT(PI»e.  1 

900  y?Y  a  WXY.-VX 

901  r  =  r*rr 

902  IF(Y.LF.Y32I  GO  To  3i 

903  31  PALL  6T0FP(i  0 , «LP • *UR • YT , 22 ) 

9O9  VJPS  a  VXY  +  Li) 

90?  RlTURXi 

906  r..p 

907  r 
906  r 

909  r 

9ln  c 

9l  1  r  ‘.UHKOOT  Ifit  uTpFP 

912  »•  9d«*0UTH  E  CALCULATES  Thf  CIV'TImuTIpN  OF  EOGF 

9l?  C  "IFF*.  ACTIO*  FPO*l  «  PLATE  hY  “CAMS  OF  Tnf 

9 1 9  c  (  tr"«tT9ltAL  ThEOMY  OF  DIFFhAfTiOf  . 

913  r  i-.u’fp  p£v.i,.i9tr.  ins  position  pf  r^t  platf  a.njo  raoak 

910  r  7.5  INPUT  ...  JIh  OUTPUT  i;t  Ihg  UO»  ThE  OlffOACTEn  fIElO. 

917  C  SUBROUTINE  VH  IS  USED  TO  CALCULATE  TML  DIFFRACTION 

9 1 3  C  FUNCTION. 

919  C 

920  r 

921  SU««OUTIr,t.  1-TnFP  (  Un  .  ALP «“UP  .  YT  ,  H  1 

922  rOM"0‘l/S"°F/PuASE  «HEI  1  .HP. 00.9.9 

92?  ru*  ?0f./A.\T/ewv.0V9.  a.*,fl.E0.?.«P 

929  <-U‘’PLt.X  CEF0,rlNC.CSPO.Vni,VP2 

925  rO“PLtt  10A.VPC.CAP 

926  TOPPLE?  J.APG.VbCl .VPC2.UU 

927  0 1 A  L  '-HI.rUK.KP 

926  J  a  (II.,  1.1 

929  tf,  a 

9  30  POP  =  YT-mF  «SI*l|i»OS  1 

931  ’itP  =  Z2»MpaCOSU*JRI 

932  P«0  =  PI  <  1.0  >*,"UH-ATAN  (  HriH/VF*  > 

933  F G  a  nrrt/SlMPl  <  1 .0  >  *91  9-PH01 

93*.  P9  a  ALP-«TA'I(V£H/mO*<1 

935  OfiOO  a  P60»?7,29b 

930  f  nTF  a  FXP < -A* I TAN( PH )) «»2 » 

937  AHrt  3  J*rF«  <‘<r./if.-2.»RGl 

93*  CAP  a  CfXP(ARG) 

939  FoH  a  90/2. /EL 

990  PrtOOO  sO. 
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•  1 '■  o  rtalb'ii  DO  »OT 

ftKPROPUCi  Li'OIiTLY. 


*♦*♦1 

44* 

i»H4 
44  3 
4  4b 

44  7 
44b 
449 
431  C 
4bl  r 

45  ?  r 
4  5  *  r 
4  5h  C 

435 

45b  r  |  m 
43  7  r ! 1  * 
43U  r  » • » 
4b*<  CM' 

4b>'  r  • 1 ! 

4bi 

4b/ 

4o  5 
4b- 
4  bs 
4  br 
*♦37 
40* 

4b** 

4  7  0 
471 
47/ 

4  7  A 
47- 

4  75 

47b 
u  7  7 
4  7  4 
479 
4hu 
461 
4  0/ 

403 

484  t  0 
4 0 r>  ?0 
4  8b 
407 
4de 

489  70 
-90 
491 
4«/ 

495 

4  94  U  0 
495  *  w 


pho02  =  c.*PHon 

CALL  (  4i>e  1  *  1  •  *fiH  tf-HOLfl  •  F'M 

r  At  L  vH  (  4*JP/  .  \|UH^«5^H,PHnLV  ,F*> ) 
»'Wri  =  C-PL^O'IHI.VJ^H 

■  :irfc*  S 

whf  =  \-ai-.nr2 

i  *u  s  4^'TT  *c  *P*V*C  /  (  /  •  *KG  ) /*  P 

4LTU«f« 

►  1*1  ’ 


cul’«i'JTlf  t  p^nELSCC *S« »S) 

T  r«  I  S  I»  Tr-t  fMtSKP.L  I'.TFr.KPl  SOo^OUTI  l  THt  l''T£GR'L  1 S  ^*0" 

•  -••*  T'J  *>,  7wr  |7'Tf««A7:0  is  r<H(  »j*P!  /<f ,  )  ,  A  ^  t*  T^t  OUTHI.T  ]S 

n»f>)«  ), 

~i«t  iston  A(i?)*Hu/)*ccd/)*r>(i/» 
ii.r.JCAI  L  Tl  ST  •LC,f-  •4' 
rg*'»41,1,  /p  i  w  P  I  •  Tk  1  •  t>'S 
r u*  *' u  *  /  T  l  s  T /lTT.ST  •  L  Jt^LG 

r  AT.\  a/1.  :  ^^7^914  o, -0. onootl  1^0^, -b,pn63t4a54t-0.noo«)7b.Ablih.  9k  0b91 

•  ^n^,«Jti  1*  ' 9**67 •- 3. 05^465040,  «o.  07‘7t>/4l  9 • 0 ./ 504037 0 1  * -0 , 1'<5b39li4 

•  i  , -0.140  c  30  9*0,1.  ♦  0  3u 40477°/ 

aT\  H/-ll#LCOOOUO33t‘*#P5?»307,;^4,-L.POUU9/9ir  ,.7. 7mU0?040  0  ,-O.OL»9b2 

•  Urtt»s,ti.L7Mfel2?b»-'M  5**m»«7t.i  .  3b?  7*  4\  ^4 ,  -  3 . 4„  *34  0/ 7b .  u  .  7n*  ^  0  \ 

•  b,-0.:M-l9‘'o/o*0.L  1^547001/ 

r  4  TA  CC/L.  ,  -  0,  0/4  9  7  3979  ,n.UP°0u3°36,".  on*  7  TT"**  ,0. Op 0689*  9?,-0.uC9 

•  4Q7 1  74  f  .  r.  ll«u 001!  9 ,  -0 . 0 0 b 74 1» 0  7  j ,  0  .  Of  n2un*2P  t  0 .  u C.  /  X  0 <  3b  7 ,  -0.00 1 2 1  7  q 

•39,0. 0*1 c3 39?o/ 

i'ATA  )/n.in«s47li4Oto.OOOnouO?A,-«.ono35i34ltO,nOl‘O^3oO6«O.RO**0^l4b 

•  ©•U.8>MVC'5c'!P#-C.*U7l??om,b*UjsnbuCb7,-O#r/7'ai0°3R,U.Olb4  975Ufi,-O 

*.'b3404M 5 *0.0009; 330b/ 

TMXS.Lt.O.C)  GO  TO  414 
VsKS 

►  z  Pl***X//.0 

ro=o.o 

C  IS 0.0 
-St  A 

1  f.  {  x  -  4  •  0  )  11,40,4  0 

v=x/4 . n 

►  SH-l 

)  )  •* 

n  =  «Pi>‘M- » ) * v 

T  p  (  *  -  /  )  >0,30*20 
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